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Phase comparators used in phase-locked loops extracting symbol timing 
from baseband data waveforms typically only produce a useful error signal 
when a data transition occurs. This gating of the error signal by data tran- 
sitions makes the natural model for studying jitter in such phase-locked 
loops time-varying and difficult to analyze. In this paper we show how, 
under good approximation, it can be simplified to a time-invariant model 
that is easily analyzed. Using this model, we study jitter accumulation 
along chains of digital repeaters with phase-locked-loop timing extrac- 
tors. A numerical example is given. 

I. INTRODUCTION 

To decode a baseband data waveform, a clock signal giving the 
proper sampling time must be available. Pilot tones are sometimes 
transmitted along with the data waveform for this purpose but, 
alternately, timing can be derived directly from the data waveform 
itself. One approach to self-timing is to let the data waveform passed 
through a memoryless nonlinearity ring a tuned circuit with a resonant 
frequency close to the nominal signaling rate. 1-8 Another technique, 
which in general involves more circuitry but gives superior per- 
formance, is to use a phase-locked loop (pll). 

The element of a pll extracting symbol timing that is most interest- 
ing functionally is its phase comparator. Numerous realizations are 
possible. Almost all these realizations are similar, however, in that 
they only produce an error signal when a data transition occurs. The 
gating of the error signal by data transitions in a pll extracting symbol 
timing makes the analysis of such a pll potentially quite different 
from the analysis of phase-locked loops used in other applications. 

For purposes such as studying timing acquisition, the effect of 
transition gating is adequately modeled as a multiplication of the 
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phase comparator gain in the presence of a data transition by the 
probability p of a data transition. After this approximation is made, 
the loop has the conventional structure. For studying jitter, however, 
this model is not satisfactory because transition gating is itself often 
an important source of jitter that of necessity is neglected when the 
average-gain approximation is made. As will be discussed more fully 
later, transition gating produces jitter by chopping the static phase 
error needed to pull the voltage-controlled oscillator (vco) from its 
natural frequency to the line frequency and thereby injecting wide- 
band noise into the loop. 

The purpose of this paper is to develop a time-invariant linear model 
for a pll with transition gating that adequately mirrors, even for jitter- 
studying purposes, the effects of transition gating. The model we obtain 
differs from the linear model obtained by simply making the average- 
gain approximation in the addition of an additive wideband noise 
representative of the chopping of the static phase error by transition 
gating. 

With the time-invariant linear model available, it is straightforward 
to calculate the spectrum and variance of the jitter on the output of 
a pll extracting symbol timing from a data waveform. We conclude 
with a discussion of jitter accumulation along chains of digital repeaters 
with pll timing extractors. 

Saltzberg 9 and Roza 10 have both previously done excellent work 
analyzing phase-locked loops extracting symbol timing. Our work is 
closer to Saltzberg's, and we use many of his results. Saltzberg was 
able to obtain equations for the jitter spectrum and jitter variance 
directly from the natural gated model. The equations we obtain from 
the approximate model are similar, but they are obtained with much 
less difficulty and without needing to assume one-sided intersymbol 
interference. 

Roza's analysis is at a more abstract level. His equations for the 
spectrum and variance of the jitter do not require an exact specification 
of the phase comparator being used as do Saltzberg's and ours but, 
consequently, they relate the jitter spectrum and variance to more 
abstract quantities with less physical meaning. 

II. GENERAL MODEL 

The general pll model we assume is shown in Fig. 1. The phase 0(0 
of the data signal and the phase 0(0 of the timing output from the 
pll are measured in slots (fractions of a symbol interval) and assumed 
to be slowly varying in comparison to a symbol interval. The loop is 
assumed to be in lock. 
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r(t)=2a k f (t-k-</)(k)) + n(t) 
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Fig. 1 — General model for a pll timing extractor. 

The input data waveform will be taken to have the form* 

r(t) = E a h f[t -k- 0(fc)] + n(t), (1) 

where n(t) is additive noise, f(t) is the standard pulse, and the {a k \ 
are the data taking on the values ±1. The normalization of the symbol 
interval to unity is without loss of generality and is equivalent to 
agreeing to measure time (as well as phase) in slots. The assumption 
of polar (±1) rather than, say, unipolar (0, 1) signaling is also without 
loss of generality. The assumption of binary data is restrictive. A 
generalization of our results for multilevel data would be nontrivial, 
but should be possible. 

Note that, in writing the output of the vco as sin [2irt — 2ird(t)~\, 
we are not precluding the possibility of the vco center frequency being 
offset from the line frequency. With the loop open and a center fre- 
quency offset of A/, the vco output is expressible as sin [2irt — 2ir6(t)~\ 
with 

d(t) = Aft + 0c (2) 

III. PHASE DETECTORS 

To proceed, we need an equation for the error signal e(t) produced 
by the phase comparator and, thus, we must specify precisely the 
phase comparator to be considered. In this section, we develop equa- 
tions for the error signals produced by two different phase comparators : 
the zero-crossing phase comparator described by Saltzberg 9 and the 
dead-zone quantizer (dzq) phase comparator used in the t4m repeater. 11 
Characterizing the dzq phase comparator will allow presenting numeri- 
cal results. 

It happens that the error signal e(t) produced by the dzq phase 
comparator is described by equations identical to those describing the 
error signal produced by the zero-crossing phase comparator if certain 
constants relating to the pulse shape f(t) and certain random variables 
relating to the channel noise n(t) are redefined. We expect that the 
error signals from other phase comparators with transition gating 
probably also have this same functional form or one very similar to it. 



* Summations and integrals written without limits are to be assumed to be from 

— oo to + °° • 



JITTER IN PHASE-LOCKED LOOPS 39 



For this reason, it is felt that, although the analysis of Sections IV 
through VII is only directly applicable for a zero-crossing phase 
comparator or a dzq phase comparator, it can probably be readily 
extended for any other phase comparator exhibiting transition gating. 

3.1 Zero-crossing detector 

The zero-crossing phase comparator produces an error signal e Z c(t) 
that is a train of weighted pulses with separation T = 1. The pulse 
weights are proportional to the time difference between zero crossings 
of the data waveform and the vco output. For a pll with a closed- 
loop response cutting off well below the symbol rate, as should always 
be the case, the exact phase of this pulse train and shape of its pulses 
are unimportant, and we can to good approximation take it to be 
given by 

ezc(t) = ZeMt-k), (3) 

where 8(t) is the Dirac delta function and the {e k \ are weights. 
A typical received data waveform, 

r(t) = E a k f(t -k-<t> k ) + n(t), (4) 

where 

4> k = <t>{k) (5) 

is drawn in Fig. 2a. It is assumed that n(t) and {<£*} are identically 
zero, and that the standard equalized pulse f(t) has a raised cosine 
shape (see Fig. 2b). Notice that, whenever there is no data transition, 
there is no zero crossing. In this case, the zero-crossing phase detector 
sets e k equal to zero. When the noise n{t) or the input phase samples 
\4> k ) are not identically zero, or when f(t) is not of its nominal raised- 
cosine shape, the received waveform is no longer as clean as that shown 
in Fig. 2a. Nonetheless, it will still maintain the same basic character 
with zero crossings occurring only with associated data transitions. 

Saltzberg cleverly shows in Ref. 9 that, to within excellent approxi- 
mation, e k is given by 

e k = aid k (6 k - <f> k - w k ). (6) 

The proportionality of e k to the phase error d k — <£* through the 
proportionality constant a x is the desired phase comparator response. 
This desired response is degraded by the on-off gating of the variable 



d k = 



1, a k t* a k+l ) 

V - (1 - a k a k+1 )/2 (7) 

0, a k = a k+ i] 



associated with data transitions and the additive disturbance w k given 
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Fig. 2 — (a) Nominal received waveform, (b) Raised-cosine pulse. 

Wk = dkVk + Clk Zl a>k-mtm, 



(8) 



where 



and 



Vk = 



*k = 



n(k + \ + <t>k) 
[/(*) - /(-*)] 

' 0, A; = 0, 1' 

f(k + |), fc ^ 0, 1 



(9) 



(10) 



[The pulse f(t) has, without loss of generating, been centered so that 
f(h) = /( — h)l- The two components a k v k and o t Enfli-m<m of u 1 * 
account for the shift in zero crossings of the received waveform by 
additive channel noise and by intersymbol interference, respectively. 
We shall not derive (6) here, since Saltzberg's derivation is quite clear 
and the derivation of 6dzq(0 m the next section and the appendix, 
being similar, presents the key ideas. 

3.2 DZQ phase detector 

The operation of the zero-crossing phase comparator is conceptually 
simple, but its implementation appears difficult. The dead-zone 
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quantizer (dzq) phase comparator, which was first proposed by W. G. 
Hammett, is more readily implemented and is being used in the T4M 
repeater. By describing and mathematically characterizing it here, we 
establish a base for numerical results to be given later. 

A block diagram of the dzq phase comparator is shown in Fig. 3a. 
Its inputs are the received waveform 



r(t) = E a k f(t -*-**) + n(t) 



(ID 



and the vco output sin \_2irt - 2x0(0]- Its only block requiring 
explanation is the dzq, which is a memoryless nonlinearity with the 
transfer function shown in Fig. 3b. 

The operation of this phase comparator is illustrated in Fig. 4, 
where waveforms at the inputs and outputs of all the blocks in Fig. 3a 
are shown. It is assumed that the data sequence is as indicated, 
n(t) = 0, <fi(t) = 0, f(t) has raised cosine shape (nominal for the T4M 
repeater), and 6(t) = £ (that is, the vco phase is lagging 45 degrees). 
The true error signal shown on the next-to-last line can be taken, for 
mathematical purposes, as the train of weighted impulses shown on the 
last line with the impulse weights equal to the integral of their associ- 
ated pulses. Notice in Fig. 4 that the phase-comparator output is 
nonzero only when there is a data transition. 

As in the zero-crossing phase comparator, the effect of additive 
channel noise and intersymbol interference is to make the weights of 
the impulses in the phase-comparator output differ from strict pro- 



r(t)=5>k f (t-M> k )+n(t) 
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Fig. 3 — (a) dzq phase comparator, (b) Instantaneous transfer function of a dead- 
zone quantizer. 
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Fig. 4 — Waveforms in dzq phase comparators. 

portionality to the difference between <f>(t) and 8(t) by a small random 
amount ^through n(t) and the data {a*}]. Defining the weights of 
the impulses in the last line of Fig. 4 as {e*}, that is, defining 



«dzq(0 = E e k 8(t — k), 
we show in the appendix that 

e k = aidk(6 k — <£t — io k ), 



(12) 



(13) 



where all symbols except the channel noise variables {v k } and the 
intersymbol interference constants {e*} are as in the discussion of the 
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zero-crossing phase comparator. The changed definitions are 



and 



«* = 1 



_ n(fe + <t» k + r~) ri(fc + 0* + t+) 
U * " (20-) ( 2 + ) 



o, ft = 0, 1 

/(* + r-)/(2r) + /(* + r^/(2/3+), fc # 0, 1 



where 



0* = /(r* - 1) - /(t±) 



(14) 



(15) 



(16) 



and r + and r~ are defined as in Fig. 5. The pulse f(t) has without 
loss of generality been centered so that 



r" + r+ = 1. 



(17) 



Thus, the only difference in the two phase comparators is in the 
precise way in which the random variables [v k ] are related to the 
additive channel noise and the constants { e k } to the pulse tails. 

IV. SWITCHED LINEAR MODEL 

Using eqs. (3) and (6) as a characterization of the phase comparator 
and realizing that the phase of the vco output is the integral of the 
signal on its input, we can model the pll of Fig. 1 as shown in Fig. 6. 
In Fig. 6, A/ is the difference between the center frequency of the 
vco and the line frequency, which we have normalized to unity. The 
function H(s) is the transfer function of the low-pass filter normalized 
so that H(0) = 1, and the open loop dc gain a is given by 

a = aia2(X3, 

where ai is the phase-comparator gain [previously defined by eq. (6)], 
a 2 is the dc gain of the low-pass filter, and a 3 is the vco gain. 

It will be convenient to mix discrete and continuous notation as in 
Fig. 6. For a precise interpretation of figures such as these, discrete 



f{t-D-f(t) 




Fig. 5 — Defining figure for t and t + . 
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Fig. 6 — Time-varying model for pll. 

labels such as 0* and Wu + <t>k are to be taken as shorthand for 
£ M(* - k) and £ (to* + 0*)5(* - k). 

The model of Fig. 6 is linear, but it varies in time as the gating 
variable d k opens and closes the loop. A direct analysis of jitter using 
this model is made difficult by the gating. 

V. TIME-INVARIANT LINEAR MODEL 

The key to obtaining a time-invariant approximation is to consider 
the nature of the signal (0* — 6 k )dk at the gate output in Fig. 6. 
Letting p denote the probability of a data transition, we have 

(0* - dk)d k = (0* - e k )v + (0* - 8k)(d k - p). (18) 

The term (0* — 9 k )p in this expansion represents the low frequency 
(in comparison to the symbol rate) response of the gate to the low 
frequency signal 0* — 6 k on its input. This signal completes the feed- 
back path. Since d k has mean p and to good approximation is not 
correlated with (0* — 6 k ), the other component is a wideband signal 
with no dc. It represents jitter-producing wideband noise injected into 
the loop by the gating. 

Conceptually, the phase error <f> k — dk has three components: a 
static offset, which we denote by n, necessary to pull the vco from its 
quiescent frequency to the line frequency ; a low-frequency component 
present when <j>(t) changes faster than the loop can track; and jitter. 
Our key approximation is 

(0*- 8 k )(dk-p) = n(dk-p). (19) 

In the pll timing extractor of the t4m repeater, y. is about 100 times 
the standard deviation of the jitter when the vco offset is a worst case 
60 parts per million. Also 0(0, which is typically jitter from previous 
repeaters, is within the closed-loop bandwidth so that it can be tracked. 
Thus, for this system approximation (19) is excellent unless the center 
frequency of the vco happens to be such that n is unusually small. 
In this case, however, approximating the wideband power injected by 
(0* — Ok) (dk — p) by the wideband power injected by 

fi(d k — p) - 
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/i(d k -p) d k w k 




Fig. 7 — Approximation to time-varying model. 

is still good since the jitter-producing noise injected into the loop by 
additive channel noise and intersymbol interference is much stronger 
than that injected by transition gating. We expect that, for most 
other well-designed phase-locked loops, approximation (19) will be 
good. 

Approximation (19) substituted in (18) gives 

(</>* - e k )d k = {4> k - 6 k )p + y.(d k - p). (20) 

With this approximation, the model of Fig. 6 becomes the time- 
invariant model of Fig. 7, which can be further simplified to the model 
of Fig. 8. 

It remains to calculate the static offset n. Refer to Fig. 6. Since 
there can be no steady-state dc signal at the integrator input, the dc 
on the output of the low-pass filter must equal -A/, and the dc on its 
input must equal — A//o. Therefore, with overscores denoting mean 
values, 



(<t> k - Qk)dk = w k d k . 



(21) 



We noted earlier that <j> k - dk and d k are approximately uncorrelated. 
Thus, 

(22) 



(</,* - d k )d k = (</>* - e k )d k = up 



and 



M = — 



A/ w k d k 



pa 



V 



(23) 



The mean w k d k is a function of the noise statistics and data statistics. 
It can be shown that if the data are independent of the noise and have 
stationary statistics, 

<Ml = 0, (24) 



4>(t) + z k 



paH(s) 



0(t) 



z k = d/p) (aj k d k +At(d k -pl) +Af/a 

Fig. 8 — Time-invariant linear model for pll. 
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in which case 

„ - - g. (25) 

pa 

The model of Fig. 8 gives the output phase 8(t) as the sum of the 
input phase 4> (t) passed through a filter with transfer function 

0( f ) = P* fl W (26) 

s + paH(s) 

and added noise 

z k = (-)[«*<** + f( d " ~ V) + ^1 (27) 

passed through the same filter. Since the model is fixed and linear, it 
is easily analyzed. Although in the remainder of the paper we concen- 
trate on finding the output-phase spectrum S$(f) for random <f>(t), 
n(t), and {a*}, exact calculations of 6(t) can be made for given <f>(t), 
n(t), and {a*}. 

VI. SPECTRUM AND VARIANCE 

The spectrum Se(f) of 6{t) can be calculated under various assump- 
tions on the statistics of <f>(t), n(t), and the data sequence {a*}. We 
make the calculation here assuming 

(i) n(t), <f>(t), and the data {a*} are statistically independent and 

stationary random processes. 
(it) The random variables {a*} are independent. 

Assumption (i) is reasonable if 4>(t) is not to represent accumulated 
jitter from previous timing extractors. 

Let p' denote the probability a data bit equals 1 and let q' = 1 — p'. 
Then 

P - 2p'q'. (28) 

Notice that the transition probability is at most \, which is achieved 
with equiprobable data. 

Assumption (i) implies that the driving processes <f>(t) and [Zk\ in 
Fig. 8 are independent.* Thus, the spectrum Se(f) of 0(0 is given by 

Se(f) = \G(f)\*S t (f) + | <?(/)!'&(/), (29) 

"The {zk\ process is a function of the {a* J and U>k) processes. Equation (9) 
(zero-crossing phase comparator) and eq. (14) (dzq phase comparator) both give 
the {vk) process as a function of both the 4>(t) and n(t) processes. Thus, the inde- 
pendence of the 0(i) and n{t) processes does not at first seem sufficient for the in- 
dependence of the <j>(t) and {vk} processes. However, the functional forms of both 
(11) and (18) are such that, in both cases, the univariate statistics of the {vk\ process 
are independent of the statistics of the <f>(l) process [the stationarity of trie n(t) 
process is used here J and moreover to within the slowly-varying <f>(t) assumption 
already made, the multivariate statistics are also. Thus, the <f>(t) and {«;*) processes 
are independent to excellent approximation. 
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where S*{f) is the spectrum of +(Q and S,(f) is the spectrum of the 
process 

z(t) = £ z fc 5(« - fc - X), (30) 

with X being an epoch-randomizing uniform [0, 1] random variable, 
the addition of which is necessary to make z(t) stationary. Assuming 
the loop to be narrowband so that the cutoff of the low-pass function 
G(f) is well within the band of frequencies over which the spectrum 
of the wide-band process z{t) is flat, we have, to excellent approximation, 

S,(f) = \G(f)\*S*(f) + |G(/)|»«Sf,(0). (31) 

Since the variance a\ of 0(t) is the integral of £<>(/), we have to this 
same approximation 

A = J \G(f) | 2 5,(/)d/ + 5.(0) J \G(f) \Hf. (32) 

The first term of (31) and the first term of (32) are associated with 
the input phase <f>(t) and will not be considered further. The more 
interesting second terms are associated with the jitter added by the 
timing extractor. Notice that S z (0) provides a measure of the timing 
extractor's performance. 

It can be shown through straightforward but tedious manipulations 

that 

aM m LCM - C(D] *fi^M - M « - M) 

p V 

+ £(€, + e_ fc - 6,_i - e*_i) a , (33) 

where C v (k) is the covariance function of the sequence {v k }. The first 
term in (33) is associated with additive channel noise and will usually 
be negligible in comparison to the remaining terms, which are associ- 
ated with the pulse shape f(t) through the intersymbol-interference 
constants { e^} and the vco offset A/ through n. 

It is straightforward to show for either phase comparator that, if 
f(t) is a symmetric pulse, then 

e k = e_,_!, k m 1, 2, • • • (34) 

and therefore 

5.(0) = [C(0) - C.(l)]/p + m 2 (2 - 3p)/p. (35) 

For unknown /x, this is certainly the best that can be done. The design 

48 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1976 



JITTER-FREE 
DIGITAL 
SOURCE 



x(t) 



REPEATER 
1 



x(t-fl,(0) 



x(t-0 N _,(t)] 



REPEATER 
N 



x(t-0 N (t)) 



Fig. 9 — Jitter accumulation along a chain of repeaters. 

conclusion to be drawn from (35) and (31), therefore, is that sym- 
metric pulse shapes are best for timing purposes. 

VII. CHAINS 

The results of the previous section are easily extended for a cascade 
of pll timing extractors. Cascades are of engineering interest for the 
study of chains of digital repeaters with pll timing extractors. In a 
cascade like that of Fig. 9, we call On(t) accumulated jitter and 



A*(0 = MO - Off-i{t) 



(36) 



alignment jitter. The name alignment jitter comes from the fact that 
A N (t) represents the amount by which the actual sampling time in 
the Nth repeater is misaligned from the ideal sampling time. 

In Fig. 9, the phase of the output of the in — l)th repeater is the 
phase of the input to the nth. This fact and the fixed linear model of 
Fig. 8 imply the model of Fig. 10, where 



Gn(s) = 



pa n H„(s) 
s + pa n H„(s) 



(37) 



This general model is fixed and linear, and thus not impossible to 
analyze by conventional methods for any parameters of interest. 
However, to derive insight and to obtain reasonably concise expres- 
sions for MO and Aat(0i ^ is necessary to make further assumptions 
that reduce the number of parameters. We assume 

(i) The additive channel noise is negligible and therefore the ran- 
dom variables {v { k n) \ may be taken as zero. 
(it) The equalized pulse shapes f„(t) at the inputs to all the re- 
peaters are identical. 
(Hi) The vco offsets A/„ are all identical. 
(iv) The filters H„(s) and open-loop dc gains a„ are all identical. 



(1/pHd k (i) k <'> + /iiil(d k -p)) + AVai 



(1/pl <d k a> k IN) + /z ,NI (d k -p» + Af N /a N 





0N<O 



Fig. 10 — Model for jitter accumulation along a chain of pll timing extractors. 
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Assumption (i) is in general reasonable. Assumptions (it) through (iv) 
are less reasonable, but are still interesting in that they are worst- 
case assumptions. The first two align the input signals in Fig. 10. 
The last gives the most chain gain at the peak frequency for a specified 
amount of peaking in each repeater. 

With assumptions (i) through (iv), the model of Fig. 10 simplifies 
to that of Fig. 11, from which it is apparent that (we are making the 
wideband noise approximation again) 



and 



Se N U) = 



St 



71=1 



&(0) 



(38) 



>*»(/) = I <?"(/) |»&(0). (39) 

Interestingly, the model of Fig. 11 is identical to that used by Byrne, 

Karafin, and Robinson 6 to analyze a chain of digital repeaters with 

tuned-circuit timing extractors. One important difference, however, 

between the two analyses is that we have available a formula for z(t) 

in terms of physically meaningful parameters. Byrne, Karafin, and 

Robinson did not have such a formula and had to experimentally 

determine S z (0). 

The variances ai* and <r 2 e „ of A N (t) and 6 N (t) are, from (38) and 

(39), given by 

<rl N = J(N)S,(0) (40) 

and 

<r 2 e K = I(N)S.(0), ( 41 ) 

where 

(42) 



J(N)= J\G»(f)\*df 



and 



I(N) = || 



n-l 



df. 



(43) 



The large N behavior of J(N) and I(N) depends critically on 
whether or not G(f) exhibits peaking (i.e., whether or not there exists 
a frequency / such that |6*(/)| > 1). When 0(f) exhibits peaking, 
it can be shown by Laplace's method (see, for example, Ref. 12 or 



z(t) = (1/p) (dKtok+Mfdk-p) ) +Af/a 





Fig. 11 — Simplified model for jitter accumulation along a chain. 
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Ref. 13) that for large N 

J (N) = J(N) = 2\G(f p )\^T/l-NL(f p )2 (44) 

and 

I(N) = 1(N) = J(N) 



QUv) 



(45) 



1 - GUp) 
where f p is the frequency at which |G'(/) | is maximum, 

L(/) = ln|G(/)|, (46) 

and L(f p ) is the second derivative of L(f) evaluated at f p . 

When there is no peaking, G(f) is everywhere less than one except 
at the origin where it must equal one. Thus, as N increases, \G N (f) | 2 
and £n=i |G n (/)| 2 become increasingly narrower in bandwidth with 
amplitudes 1 and JV 2 respectively at zero frequency. Let 

(?(/) - exp{L(/) + j0(/)} (47) 

with L(/) and <^(/) real, and approximate L(f) and <j>(f) near the 
origin by 

L(f) = L(0)+L(0)/+ (h)L(0)f 

= ti)L(0)P (48) 



and 



4(f) = 0(0) + 4(0)/ 

= 4(0)/. (49) 



The quantities L(0) and L(0) equal zero because G(f) equals one at 
the origin and \G(f)\ is even. The phase 0(0) equals zero since (7(0) 
is real. Equations (48) and (49) lead after much manipulation to the 
large N approximations 

J(N) = >/*./#[- L(0)] (50) 

and 

I(N) = 2irN/<i>(0). (51) 

If (?(/) is a single-pole characteristic, that is, if 

G(/) = r+W (52) 

then (68) reduces to 

I(N) = 2rNf , (53) 

which is consistent with a result in Ref. 6. 

Notice the great difference in the asymptotic behavior of J(N) and 
I(N) in the two situations. With peaking, J(N) and I(N) are asymp- 
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totically exponential [eqs. (44) and (45)], whereas without peaking 
J(N) actually decreases with N for large N and I(N) only grows 
linearly. The rapid growth of (44) and (45) with N indicates that, to 
control jitter accumulation along a long chain, there must be very 
little if any peaking. 

VIII. NUMERICAL RESULTS 

We have used the theory developed here to numerically characterize 
the jitter performance of the timing extractor in the t4m digital re- 
peater. The exact numbers obtained are, of course, only of direct 
interest to designers of this repeater, but their relative sizes should 
be indicative of other timing extractors, and thus of more general 
interest. 

Data on the t4m line are scrambled so that an assumption of inde- 
pendent data is justified, and, furthermore, the probability p of a 
data transition equals \. Additive channel noise is a negligible jitter 
source. The timing loop has a worst-case vco offset of 60 parts per 
million, a dc gain 

a = 8 X 10- 3 , (54) 

and a filter 

„ (n _ a + jf/mi + if/f*) (m 

U) " (i + j///.)U + i/// 4 ) 2 d + of /Mi + ////•)' K } 



where/i,/ 2 ,/ 3 ,/4,/5, and/ 6 are, respectively, 3.65 X lO" 6 , 3.65 X 10" 5 , 
1.82 X 10- 7 , 7.30 X 10- 5 , 1.45 X 10~ 4 , and 1.09 X 10~ 3 cycles per 
slot.' 

The reasons for the rather involved form of this filter are of no 
consequence here. However, it should be pointed out that some of the 
poles and zeros are fixed by circuit parasitics and are not introduced 
intentionally. 

The shape /(/) of the received and equalized standard pulse depends 
on repeater spacing and other parameters. A typical waveform f(t) 
appears in Fig. 12. With a slicing level A fixed at f the maximum of 
1/(0 - f(i — l)l> tne associated intersymbol interference constants 
{ e*} are as in Table I. Substituting these values and p = \ in (33) gives 

5,(0) = n 2 - 2^(0.0090) + 1.61 X 10- 4 . (56) 



For | A/ 1 ^ 60 X lO" 6 , we have 



M = 



pa 



^ 0.015. (57) 



* Assuming T = 1 is equivalent to measuring time in slots and thus frequency 
in cycles per slot. 
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ARBITRARY TIME ORIGIN 

AMPLITUDE SCALED SO THAT PEAK IS 1 




.£ 



-0.5 - 



t (SLOTS) 



Fig. 12 — Output of differentiator. 

With m £ [-0.015, 0.015], 

&(0) G [80 X 10- 6 , 656 X 10" 6 ]. 

The minimum is achieved by 

v = +0.009 

and the maximum by 

u = -0.015. 



(58) 



(59) 



(60) 



All three terms of (33) are of comparable magnitude for these parame- 
ters, indicating that for this particular repeater the jitter produced 
by intersymbol interference is comparable to that produced by tran- 
sition gating. 

Table I — Numerical values for the constants {« fc } having 

magnitudes greater than 0.00005. The constants 

to and t_i equal zero by definition 



«-l 


-0.0009 


>■ I 


0.0194 


e-i 


0.0000 


«• 


0.0000 


<i 


0.0284 


<J 


-0.0002 


C| 


0.0003 


•4 


-0.0005 


«6 


0.0001 
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Table II — 


Integral values and their 


approximations in dB 


N 


f(N) 


J(N) 


f(N) 


nm 


1 


-40.0 


-38.2 


-27.2 


-38.2 


2 


-40.7 


-38.9 


-28.0 


-34.0 


3 


-40.9 


-39.0 


-28.1 


-31.6 


5 


-40.6 


-38.8 


-27.8 


-28.6 


7 


-39.8 


-38.3 


-27.1 


-26.4 


10 


-38.4 


-37.2 


-25.7 


-23.8 


20 


-32.7 


-32.1 


-19.9 


-17.2 


30 


-26.3 


-26.0 


-13.6 


-11.4 


50 


-12.9 


-12.7 


-0.2 


+0.9 


70 


0.8 


1.0 


13.6 


14.2 


100 


21.8 


21.9 


34.6 


35.0 


200 


92.8 


92.9 


105.6 


105.8 


300 


164.4 


164.5 


177.2 


177.3 


500 


308.4 


308.4 


321.1 


321.2 



With p = i a = 0.008, and H(f) as given by (55), (?(/) has 0.75 dB 
of peaking. In Table II we give J(N), J(N), t(N), I(N), a\s, and 
a 2 , K [as denned by eqs. (44), (42), (45), (43), (40), and (41)] for various 
values of N. The entries for a\ N and o 2 6n assume a worst-case 5,(0) of 
6.56 X 10~ 4 . The convergence of J{N) to J(N) and f(N) to I(N) 
for large N is readily apparent. It is also apparent that jitter accumu- 
lation limits line lengths on this particular system to about 70 re- 
peaters where the rms alignment jitter is 0.03 slot (10 degrees). The 
entries in the table for larger N need qualification. One could not 
expect to measure them since, in that length line, lock would be lost. 
Our analysis assumes lock and makes many small signal approximations. 

IX. CONCLUSIONS 

By chopping the phase error <j> k — 8 k , transition gating injects 
jitter-producing wideband noise into a pll extracting symbol timing 
from a baseband data waveform. If the wideband noise injected by 
the signal (<f> k - k )(d k — p) is approximated as that injected by 
n(d k — p), the model for analyzing jitter in a pll with transition gating 
becomes time invariant and can be easily analyzed. 
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APPENDIX 

Calculation of Error Signal 

In this appendix, we find the weights {e k } for the dzq phase com- 
parator as a function of the input phase, the vco phase, the data se- 
quence {a,k\, the standard pulse shape f(t), and the additive channel 
noise n(t). The procedure followed is quite similar to that used by 
Saltzberg 9 to find the weights [e k \ for the zero-crossing phase 
comparator. 

Assume first, for argument's sake, that a* = — 1 and a k +i = + 1. 
Then 

r(t) = -f(t -k- 0*) + /(*-*-!- 0*+i) 

+ £ a m f(t -m- W ) + n(t). (61) 

The input jitter 4> k contains only low frequencies in comparison to 
the symbol rate. Thus, the near neighbors of <t> k approximately equal 
it and, assuming f(t) is insignificantly small far out on its tails (as it 
must be for the eye to be open and regeneration possible), we can 
reasonably take r(t) for t G [&, k + 1] as 

r(t) = -fit - k - <f> k ) + f(t - k - 1 - <t> k ) 

+ £ a m f(t - m - **) + n(t). (62) 



m^k, k +l 



Define 



Then 



q(t) = L aj{t - m - <f> k ) + n(t). (63) 

m?t k , k +l 



r{t) = -fit- k - **) + f{t - k - 1 - 0*) + qit). (64) 

Let f denote the time in [0, 1] when fit — 1) - f(t) makes an 
upward crossing of the dzq slicing level A and let t + denote the time 
in [0, 1] when it makes a downward crossing (see Fig. 5). During 
the time interval [k + **, k + 1 + <£*], r(t) will appear as in Fig. 13. 
Define r k so that the input to the dead-zone quantizer makes an up- 
ward crossing of the slicing level A (that is, a pulse starts at the dzq 
output) at time A; + <f> k + rf . Denote the slope of f{t — 1) — fit) 
at time r~ by 0~. That is, let 

fi- = f^- - 1) - f(r-). (65) 
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Ht)-q(t) 



q(k + k +T~)| 



tq(k + k +T + ) 




4 \ 4 4 _ 

,' S k+0 k +T" / 'k+0 k +T+ 

k + k +T k - k + k +T + 

Fig. 13 — Linear extrapolation for r~ and t + . 



To good approximation, we have (see Fig. 13) 



. q(k + </>* + T~) 

r = r- - r." 



(66) 



or, after rearrangement, 

rf = r- - g(fc + 0* + r-)//3~. (67) 

Denning r^" as the time when the pulse at the output of the dead-zone 
quantizer ends and denning 

0+ = /( T + - l) - /( T +), (68) 

we have, by a similar argument, 

ft = r+ - q(k + ft + r+)/0+. (69) 

The output of the vco is sin [2?r* - 2x0(0]. Assuming the phase- 
locked loop is in lock, that is, that 0(0 - <f>(t) is, at most, a small 
fraction of a slot, a downward crossing of zero by sin [2irt — 2ir0(0] 
must occur near the middle of the interval [k + <f>k, k + 1 -\- 0*]. 
Define Sk so that 

k + <f> k + § + s k 

denotes the time of this occurrence. Then 

[2ir(A + 0* + i + «*) - 2tt0(/c + 0* + i + s*)] mod 2tt = r. (70) 

The phase 0(0 is slowly varying with respect to the symbol interval. 
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Therefore, 

6(k + 4> k + \ + s k ) = d(k) = d k . (71) 

With this approximation, (70) becomes 

[2tt(A; + ** + * + s*) - 2tt0(A;)] mod 2* = r, (72) 

which implies (4> k — 6 k and s* are small) 

s k = 6 k - fc . (73) 

From Fig. (14), it is clear that e k must be proportional to 

* + 8 k -rr- Lrt - (I + »*)], (74) 

which simplifies to 

1 + 2s* - u ~ rt (75) 

or, after substituting (74) for s k , 

2(0* -e k )-\-l- r t ~ - rt. (76) 

Using the approximations (67) and (69) for ri" and rif" and denoting 
the constant of proportionality relating e k and (76) by ai/2, we have 



+ *i\ qik+ fL + ^ ) + qik +£ +t¥) \- (77) 



e k = ai(6 k - 0*) + «i(l - r~ - t+)/2 

2/3~ ' 2tf 

Let 

"* = 2T + W (?8) 

and 

_ f(k + r-) /(fc + r+) . 

ek ~ 2/3" + 2/3+ (79) 

and assume that the pulse /(<) has without loss of generality been 
centered so that 

T - + T + = 1. (80) 



/ 


TPUT OF IDEAL LIMITER ^ — C 
, t 


UTPUT OF RECTIFIER 




1 


r n 
i 

i 

i 







k + k k + </> k +T k _ k + «/> k + 1/2 + s k k + # k +T k + k+0 k + 1 

Fig. 14 — Inputs to multiplier. 
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Then, using (63), we obtain 

e k = ai(6 k — <f> k + v k + Z a m e k - m ). (81) 

m^ k , k +l 

To derive (81), we assumed a k equal to —1 and a k +i equal to +1. 
If the signs are reversed, a similar analysis shows 



e k = ai(6 k — <f> k — v k — Y. a m e k - m ). 

m?± k , k +l 


(82) 


If a k = a k+ i, 

e k = 0, 


(83) 


since no transition occurs. By defining 




d f 1 ' *"* ak+1 ] _ (1 - a k a k+1 ) 
10, a k = a k+ ii 


(84) 



we can combine all three of these equations, and write for all a k , a k +i 
e k = aid k (8 k — <f> k — a k v k — a k £ a m e*_ m ). (85) 
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